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^vq We abstract the essential features of holographic dimer models, and develop sev- 

en 

eral new applications of these models. First, semi-holographically coupling free band 
fermions to holographic dimers, we uncover novel phase transitions between conven- 
tional Fermi liquids and non-Fermi liquids, accompanied by a change in the structure 
of the Fermi surface. Second, we make dimer vibrations propagate through the whole 
crystal by way of double trace deformations, obtaining nontrivial band structure. In 
a simple toy model, the topology of the band structure experiences an interesting 
reorganization as we vary the strength of the double trace deformations. Finally, we 
develop tools that would allow one to build, in a bottom-up fashion, a holographic 
avatar of the Hubbard model. 
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I. INTRODUCTION 

Holographic models of condensed matter systems have seen a recent surge in interest. 
While this program can by now point to a few successes, there is one feature of realistic 
solids that is commonly not shared by their holographic stand-ins: in real solids, translation 
invariance is broken to a discrete subgroup by formation of a lattice. 1 In the holographic 
models with translational symmetry unbroken, momentum is strictly conserved as there 
are no Umklapp processes to dissipate it. Consequently, in the background of an electric 
field, energy and (at finite density) momentum are pumped into the system at a constant 
rate. This leads to unrealistic transport properties such as an unsmeared delta-function 
Drude peak, among other things. 2 Also, in heavy fermion materials, strong correlations 
of itinerant electrons with localized spins on a lattice are believed to trigger interesting 
phenomena such as quantum criticality. It is desirable to find holographic models which 
exhibit lattice structure. 

In [BJ, we constructed holographic models with translational symmetry explicitly bro- 
ken by fermionic degrees of freedom localized on lattice sites, interacting with a continuum 
gauge field. 3 These systems have many appealing features. They naturally give rise to bulk 
fermions living on a lattice of AdS2 spacetimes; such spacetimes play a crucial role in holo- 
graphic non- Fermi liquids [T014T3] . Thus they provide a natural home to study holographic 
non-Fermi liquids without worrying about the physics of the asymptotically anti-de Sitter 
(AdS) Reissner-Nordstrom black brane, especially its large ground state degeneracy and 
potential instabilities. It has also been argued from the field theory point of view [2] that 
a lattice of localized defect fermions interacting with continuum degrees of freedom is the 
best candidate for a condensed matter system exhibiting the phenomenology of holographic 

1 A second, somewhat related, feature which is lacking in most holographic studies of condensed matter 
systems is disorder. For early studies of how to include disorder, see [HE]- 

2 One way to work around this is to make use of the large number of internal degrees of freedom, which is 
one of the defining features of holographic models. A plasma sea with order A'' 2 excitations, where A is a 
large integer, can effectively act as a heat and momentum dump for charge carriers, leading to interesting 
DC conductivities [3H5]- 

3 An alternative proposal in which the whole gauge theory was forced to live on a lattice was put forward in 
[7J. A method for incorporating lattice defect fermions through semi-holographic techniques, which works 
only when the defect fermions are neutral under the large N gauge group, appears in [5]. Other lattice 
defect models with related features to the models we study will appear in [9]. 



3 



non-Fermi liquids 

In this paper, we abstract a few essential features of the holographic dimer models con- 
structed in [6], and then extend our knowledge of these models in several directions. First, 
we add free band fermions to the boundary theory and then weakly mix them with the 
large N sector semi-holographically, following Faulkner and Polchinski [8]. When the large 
N sector undergoes a dimerization transition, it induces a dramatic change in the singlet 
fermion sector. Namely, the melting of dimers turns a normal Fermi liquid into a non- Fermi 
liquid, with an accompanying change in the structure of its Fermi surface. This transition 
may be of some interest in relation to toy models of heavy fermion materials [H] . 

We then turn to a study of dimer vibrations. In the original holographic dimer models, 
the dimers do not effectively talk to each other in the large N limit, and each has a discrete 
vibration spectrum. By adding double trace deformations, we let the dimers communicate 
with each other, allowing their vibrations to propagate through the whole crystal in the form 
of Bloch waves. As an illustration, we explicitly work out the band structure in a solvable 
toy model of this type. Interestingly, we observe a reorganization in the topology of the 
band structure as we vary the strength of the double trace deformations. 

Finally, we also outline how one would deform the system by an actual fermion hopping 
term so as to liberate otherwise immobile localized defect fermions. The corresponding 
problem in the bulk is beyond current brane technology, as one would have to know the 
non-Abelian generalization of the Dirac-Born-Infeld (DBI) action describing the fluctuations 
of the brane. Only the lowest few terms in a derivative expansion of this action are currently 
understood. However the deformation we identify can easily be implemented in a bottom- 
up model, where the bulk theory is not taken to be a system realized in string theory 
but instead is governed by an effective action with free parameters which can be matched 
against "experiments." With this toolkit it is in principle possible to construct holographic 
avatars of one's favorite lattice models such as the Hubbard model. This should allow one to 
parametrize his/her ignorance in terms of the coefficients in the bulk effective action. While 
such a program has been somewhat successful for QCD [151 IB], we will refrain from any 
such attempts in this work. 

The organization of this paper is as follows. We first abstract out the salient aspects of 



our previous construction [B] in Sec. [TTJ In Sec. Ill we then describe semi-holographic phase 
transitions between Fermi liquids and non- Fermi liquids, much in the spirit of [5]. Changing 
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FIG. 1: (a) The high-temperature phase with U(l) x U(l) symmetry. (b)The low-temperature 
phase with {7(1) x U(l) —> {7(1) via brane recombination. 



gears, in Sec. |IV[ we show how one can map out the band structure of dimer vibrations in 
the holographic dimer models perturbed by simple double trace deformations. For a soluble 
toy model, we observe reorganization in the topology of the band structure as the strength of 
the double trace deformations increases. Lastly, in Sec. [Vj we outline how knowledge of the 
non-Abelian DBI action would enable us to study a genuine fermion hopping deformation 
and how one could go about building holographic bottom-up models of generic strongly 



coupled lattices. Detailed calculations supporting the plots of band structure in Sec. [IV 
have been relegated to the appendix, "Hypergeometric-ology." 



II. PHYSICS OF A HOLOGRAPHIC DIMER 



We first review the essential picture of a single dimer in the holographic dimer lattice mod- 
els constructed in [B]. We began with a pair consisting of a D5-brane and an anti-D5-brane, 
with the bulk spacetime being the AdS5 Schwarzschild black brane [see Fig. 1(a)]. Asymp- 
totically, they wrap copies of AdS2 x S 4 in the AdSs x S 5 . In the boundary field theory, they 
introduce localized defect fermions, transforming under fundamental and antifundamental 
representations of the SU(N) gauge group of the continuum gauge field, respectively. As 
we cool down the system, the size of the black brane shrinks down, and the D5-brane and 
anti-D5-brane pair up [see Fig.l(b)], spontaneously breaking {7(1) x {7(1) down to a diago- 
nal i7(l). In the boundary field theory, this transition is characterized by a mildly nonlocal 
order parameter involving fermions from neighboring sites as well as an open Wilson line to 
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insure gauge invariance [T7]. This is the essence of the dimerization transition worked out 
in detail in [UJ. 

We now look, rather abstractly, at fluctuations of these probe branes. 

A. Undimerized phase: gapless spectrum 

In this phase, we can study each probe brane separately. For every possible fluctuation 
of the probe brane, there exists a corresponding gauge singlet operator localized at a point, 
Oj, bosonic or fermionic. Here J labels whether the operator is associated with a D5-brane 
or anti-D5-brane, and it will be promoted to a lattice index when we discuss a lattice of 
dimers. 

By studying fluctuations in the bulk, it is in principle possible to work out 

J dte^(Oj(t)O\,(0)) = i6 Jt j,g{u>) (1) 

in detail. For example, had we stayed in the undimerized phase down to zero temperature, 
the induced metric on the (anti-)D5-brane is exactly that of AdS2 and (0 + l)-dimensional 
conformality would dictate 

G{u) = ecu 2 *- 1 (2) 

with c a calculable complex number and A the operator dimension of Oj. 4 

When we turn to semi-holographic constructions, it is crucial to keep in mind that Q{oS) 
in our concrete holographic lattice model behaves differently from the Green's function for 
AdS2 x R d_1 at finite temperature, appearing in the construction of holographic non-Fermi 
liquids. This is because the metric on the embedded D5-brane, induced from the AdSs 
Schwarzschild black brane metric, is different from the near-horizon geometry of the AdS 
Reissner-Nordstrom black brane. 

4 Strictly speaking, we can only trust our analysis in the undimerized phase down to temperatures of order 
T 3 a 3 ~ y/X/N. At this point the backreaction of the D5-branes on the background geometry can no longer 
be neglected [BJ . For most of this work we can safely neglect this complication as we are shielded from these 
parametrically small temperatures by the dimerization phase transition which occurs at Ta ~ 1. However, 
if we considered a lattice made purely of D5-branes instead of alternating D5- and anti-D5-branes (as was 
also discussed in [BJ), there would be no dimerization transition. In that case the physics would be well 
captured by the AdS2 gravity background down to this parametrically very low temperature. 
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Nevertheless, one generic behavior of this type of model is that, since the probe brane is 
touching the horizon in this phase, the spectrum is gapless. This means that lim w ^ Gif^) — 
0. 



B. Dimerized phase: gapped spectrum 

In this phase, probe branes recombine and no longer stretch down to the horizon, leading 
to a gapped spectrum. Still, the induced metric near the asymptotic boundary is that of 
AdS2- Let us write the asymptotic AdS2 metric as 

rf 4ds 2 = i(~ d * 2 + dz 2 j). (3) 

z j 

Then near the boundary point where the anti-D5-brane stretches down (zj^ = 0), a time- 
dependent fluctuation with frequency co can be expanded as 



g-iwt 



«d» {^ A + -} + ftsK") i z m + -}] ( 4 ) 



whereas near the boundary point where the connected D5-brane stretches up (zds = 0), the 
same fluctuation can be expanded as 

e-** [a m (u) {z^ A + ...}+ fa( u ) {z£ 5 + •••}] • ( 5 ) 

The coefficients a^(co), 0f^(cu), «D5(w), and /5ds(w) are related by a frequency- dependent 
matrix as follows: 

hi(u) t 22 (u) ) \/3 m (u) ) ~ \Pm(e>) ) ' 

For a simple phenomenological toy model, T(u) is worked out in detail in Appendix [Aj 

In the absence of any deformation of the theory, there exists a steady dimer vibration 
with frequency u n if and only if there exists a consistent nontrivial solution with a^^n) = 
«D5(^n) = 0. This can happen if and only if t 12 (u n ) = 0. This gives rise to a pole in 



/ a D5 (uj) 1 









fdte iu]t (Oj(t)O],(0)) at 

Typically there exists a gap to the first dimer excitation and at low frequency correlation 
functions give 

lim J dte^(Oj(t)O^(0)) = iA JtJ ,, (7) 
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with all the components of Ajji generically nonzero. 

In passing, we note that for the top-down D3/D5 system of [6], the worldvolume gauge 
field and the slipping mode scalar mix due to the Wess-Zumino terms in the action. The 
mixed sector gives rise to two towers of scalar fields depending on the angular momentum I 
on the internal sphere. As shown in [18] they behave like fields with mf = (I + 3) (7 + 4) and 
mf = 1(1 — 1), respectively. At / = we are hence effectively describing a massless scalar, 
presumably dual to the defect fermion bilinear, as well as a massive scalar mode dual to a 
dimension 4 operator. 

Now that we have gathered essential information regarding the dimerization transition 
and spectra in both phases, let us look at several physical applications. Anticipating the huge 
landscape of large N dimer models, we will keep A and all the other information (computable 
in explicit models) as undetermined free parameters. From here on, each subsequent section 
of the paper can be read independently. 

III. SEMI-HOLOGRAPHIC PHASE TRANSITIONS 

A rich set of non-Fermi liquid behaviors has recently been discovered by studying the 
physics of probe fermions in the asymptotically AdS Reissner-Nordstrom background [TOl - 
[13] . The near-horizon AdS2 X R 2 region of the black brane plays a crucial role in organizing 
and explaining this physics; the physics of the emergent "locally quantum critical" theory 
dual to the AdS2 region is what gives rise to the non-Fermi liquid behavior. However, this 
black brane is in some ways nongeneric. For instance, it suffers from a superconducting 
instability in the presence of generic charged scalar fields in the bulk [TJJ], and neutral 
scalar fields coupled to the kinetic term of the bulk U(l) gauge field deform the near- 
horizon geometry [20] to be of the Lifshitz form |21| . Even the backreaction of the fermions 
themselves deforms the near-horizon geometry to Lifshitz form at sub leading orders in 1/N 
[22] [shifting the AdS2 x R 2 geometry, which has a dynamical critical exponent z = oo, 
to instead have z ~ N]. While in many cases these deformations may leave the essential 
physics of the fermion spectral function unchanged (see [8 J for a nice discussion), it is also 
reasonable to find other ways that the essential insights of [Tom3] can be reproduced in a 
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more robust setting. The A0IS2 regions spanned by the D5- and anti-D5-branes in the top- 
down holographic dimer model of [H] provide an alternative way to obtain the same physics. 
Here, we explore this in a semi- holographic setting following [8|, and we abstract the main 
features of the top-down model to include more generic possibilities. 

We begin with a large N field theory, governed by some action Strong, with the following 
features: 

1. There is a lattice of defect fermions which undergoes a dimerization transition as we 
vary the external parameters such as temperature [see Fig. 2 for the (l+l)-dimensional 
case]. We will focus on the cases for which this parameter is temperature, but one can 
easily generalize. 5 

2. There exist fermionic operators Oj localized at the Jth lattice site, whose thermal 
correlation functions in the undimerized phase are known and gapless: 

J dte^(O F (t)O^(0)) = i5 JtJ ,g(u), 

with g(u) ~ w 2A_1 for u > T. (8) 

3. In the dimerized phase, the spectrum is gapped and 

lim J dte iujt (O F (t)O F }(0)) = iAj t j>. (9) 

Here, Aj t j/ is nonzero (generically if and) only if J = J' or J and J' are paired up via 
dimerization. 

For example, for the literal D5 probe theory in AdSs x S 5 , we can take 

o F = xW=Mxj (10) 

and work out G(oS) and Aj t j> as a function of external parameters. Here Xj^ = i(J) is the 
Af = 4 gaugino evaluated at the Jth lattice site, and \j is t ne probe fermion associated 

5 For instance, one can consider driving such a transition by going to finite chemical potential for the large 
N gauge fields at T = 0, at the cost of introducing Reissner-Nordstrom black branes. At sufficiently 
large chemical potential, even at zero temperature, the horizon of the extremal Reissner-Nordstrom black 
brane grows large, and the probe branes will transition back to a configuration where they stretch to the 
horizon instead of reconnecting. It would be interesting to determine the order of this phase transition at 
zero temperature. 
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with the Jth site. There is also an infinite tower of similar operators of higher conformal 
dimension. We will, however, keep our discussion abstract. 

Note that we are only guaranteed of the scaling form ^ governed by the (0+1)- 
dimensional conformal invariance when u ^> T. 6 Therefore, looking forward for a moment 
(to the stage where we mix the O f s with semi-holographic fermions) this behavior of the 
Green's function will be relevant when studying excitations close to the Fermi surface, only 
if the disconnected phase persists to very low temperatures (compared to the Fermi mo- 
mentum kp). This is achievable in our models, because the temperature of the dimerization 
transition is T c ~ 1 ^ [B] , where affect is the lattice spacing for defect fermions, and can be 
dialed freely; while kp ~ — - — is another free parameter, where aitmerant is the lattice spac- 

^itinerant 

ing for semi-holographic itinerant free fermions, which can also be adjusted independently. 
Thus we make a hierarchy affect 3> ^itinerant- 



We now semi-holographically couple this large N field theory to the free band fermion in 
the spirit of [S]: 7 

S = Strong + ^ / dt Cj( iS J,J' d t + ^J, J' + tj,j')cj' 

J J 

Here tj t j> characterizes the band structure of the originally free fermion c sector, which now 
mixes with the large N dimer model through the coupling constant g. 

The key insight of [8] is that large N factorization of the field theory (which would work 
even at small 't Hooft coupling) can be used to infer the modifications to the two-point 
functions of the conducting c fermions arising from the coupling g. The g = Green's 
function for the c fermions is 

Go(k..) . -iJ- E /, ie .-«^> (cj(i)c », ( 0)) 9rf . _ i , |k 1 _ kf(k)| (12) 



6 At very low frequency G(^) will still approach zero on general grounds, but it may do so with a different 
scaling dimension A' or in even more complicated ways. 

7 For notational simplicity, we made the free c fermions live on the same lattice sites as the defect fermions 
do. As just mentioned, however, we should really make the c fermions live on a much finer lattice to get 
the hierarchy — ~ kp ^> T c ~ — . Also as in §2 of \8[, we have neglected possible spin-orbit 



effects that could promote, for example, coupling constants between c and O f in (11) to be matrices in 
spin space. 
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with kp(k) the point on the Fermi surface, closest to the argument k, and N\_ 8 , the number 
of lattice sites. Then we find that for finite coupling g, after summing a geometric series of 
tree-level mixing diagrams, 

G g (k,u) ,\. ^pq- r . (13) 

to — v\k — k F {k)\ — g 2 y{i£, u) 

In particular, for Q{k, u) = <%o> 2A_1 with A < 1, one finds a dominant low- frequency behavior 
characteristic of a non-Fermi liquid which has vanishing quasiparticle residue [with marginal 
Fermi liquid behavior precisely at A = 1, when the naive a; 2-1 is modified to have a;log(a;) 
behaviour]. For A > 1, the residue does not vanish, but the theory is still novel in that 
the quasiparticle width does not agree with that of standard Fermi liquid theory. As we 
described above, these results are true in a regime where kp 3> u ^> — - — , where the 

^defect 

zero-temperature Green's functions used above should be a good approximation to the true 
(finite- but low-temperature) answers. 

Now, we are in a position to add one simple observation on top of the basic picture 
advocated in [8]: in holographic models which undergo a dimerization transition as in 
Sec. [IIJ the phase transition also drives an interesting transition in the structure of the 
Fermi surface. The main point is that the low frequency behavior of the Green's function 
J dte lU)t (Oj (t)Oj) (0)) changes drastically in the dimerization transition. In the undimerized 
state, we will have non- Fermi liquid behavior just as in [Sj. However, in the dimerized phase, 
the spectrum in the dimer sector is gapped. This means that at low frequencies, instead 
of exhibiting power-law behavior, lim^o £?(k, ui) = A for some nonzero constant A. Thus 
in this phase, we have a conventional Fermi liquid whose Fermi surface is shifted from the 
original kp. 

Therefore, in this semi-holographic setting, the dimerization transition of Sec. [TT| becomes 
a transition between a conventional Fermi liquid phase (dimerized) and a non-Fermi liquid 
phase (undimerized). These transitions are somewhat reminiscent of the phase transitions 
in Kondo lattice models discussed in [14] and references therein. 

Finally, we note that if one is purely interested in finding realizations of the non-Fermi 
liquid phase, without studying phase transitions of the Fermi surface, one can also simply 
study the BPS lattice model made only of D5-branes and generalizations thereof. In this 
case, the Schwarzschild AdS^ black brane (with probe D5-branes wrapping AdS2 subspaces 
in the AdSs) correctly captures the physics down to temperatures of order T ~ n 1 ^*' — 



^defect 
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FIG. 2: Dimerized configuration of interest. Note that we made the distance between J = (2j+l)th 
site and J = (2j + 2)th site smaller than that between J = (2j)th site and J = (2j + l)th site, so 
that we have one unique dimerized configuration below the critical temperature. 

After incorporating semi-holographic fermions, now without the constraint uj 3> — - — as 

^defect 

there is no dimerization transition, even very low frequency excitations above the Fermi 



surface are governed by (\8\) and (13) at large N. 



IV. DOUBLE TRACE DEFORMATION: BAND STRUCTURE OF DIMER 

VIBRATIONS 

Next, we visit the landscape of holographic dimer models with certain double trace defor- 
mations added to the boundary Lagrangian. For simplicity, let us consider the 1-dimensional 
array of dimers (see Fig.2). Let us label sites so that the J = (2j + l)th site is paired up 
with the J = (2j + 2)th site with j 6 Z. Our inputs are: 

1. There are bosonic Hermitian operators Oj which corresponds to a bosonic fluctuation 
of a probe brane originating from the Jth site. The fluctuation take the asymptotic 
form (for frequency uj) 

e-** [aj{u) {zy A + ...}+ {4 + ...}] • (14) 

2. We stay in the dimerized phase where the coefficients a 2 j+i(w), /3 2 j+i(u), a 2 j+2(^), 
and /%+2(w) are related as follows: 

«2j+2(w) \ _ I t 12 {u) \ I a 2j+ i(u) 

/?2i+2(w) / I *2l(w) t 22 (uj) J \ 
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Originally, the dimers are basically decoupled from each other and each has a discrete 
vibration spectrum at w = w„ where £12(^1) = 0. We now deform the theory by double 
trace operators, and determine the resulting band structure. 

A. Double trace deformation 

We add a double trace deformation of the form 

AL A ,. = h'J20*0* +1 (16) 

to the Lagrangian. The effect of double trace deformations on the dual gravitational de- 
scription is well known [23H27] . In our context, the standard recipe leads to 




(17) 



where h = (2A — l)h'. Note that this is a relation between (2j)th site and (2j + l)th site 
belonging to different dimers, not a relation between the (2j + l)th site and the (2j + 2)th 
site which together form a dimer. 

Instead of the totally reflecting boundary conditions implied by the undeformed condition 
aj = 0, the deformed boundary condition ( JT7[ ) implies that when a fluctuation associated 
with (the bulk field dual to) hits the (2j)th AdS2 boundary, part of the wave gets 
reflected back towards the (2j — l)th site, but a small fraction (governed by h) instead gets 
transmitted to the (2j + l)th site. 8 

B. Band structure 

In summary, the equation of motion on the probe brane (top-down or bottom-up) gives 



the relation (15) whereas the double trace deformation yields the relation (17). We now look 
for most general (spatially normalizable) time-dependent, but nondissipative, modes with 
these relations. For simplicity, we use the periodic boundary condition with a number of 
lattice sites N\ St = 2N dimer and then take the thermodynamic limit N dimer — ?• 00 at the end. 



These "transparent" boundary conditions modify the propagator of the scalar fields with interesting 
consequences for loop corrections, as explored in 
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Starting with a generic (ai, (3i) and evolving through the chain, we get: 

T(w) 







V*) 













*2l(w) *22(w) 



)l 










/i 



i / V *2l(w) *22(W) 




(18) 
(19) 



Continuing this way, and getting back to the original site, our periodic boundary condition 
imposes 




ht 2 i(uj) ht 22 (ui) 



h 



h 



)] 


-^dimer 

1 











(20) 



Taking the thermodynamic limit iVdi mer — > oo, we conclude that there exists a nondissipative 
solution with frequency u(k) and with Bloch momentum k if and only if 



det 



ht 2 i{uj{k)) ht 22 {u(k)) 

tuMkft t 12 (io(k)) 







Aka 



h h 

The band structure is encoded in u(k). 



with k G 



7T 71 



a a 



(21) 



C. A simple toy model 

If we know T(co>), it is a simple matter to map out the band structure numerically. As 
an illustration, let us perform this exercise for the caricature toy model described fully in 
Appendix |A} For the special case of A = 1 we display the first few bands for several values 
of h in Fig. |3j 9 Without the double trace deformation, the full spectrum is just N^ mpr 
copies of the spectrum of a single dimer, independent of k. We see that for small h we still 
have a band structure with very narrow bands centered around the mode spectrum of the 
uncoupled dimers, u n = — for nonzero integer n. Around h = 0.8 a new band emerges and 
this sector becomes gapless, potentially signaling the onset of an instability towards forming 
a spatially inhomogeneous condensate similar to the one encountered at the Lifshitz point in 
ferromagnets. As h is the coupling constant of a double trace deformation, such a spatially 
inhomogeneous condensate would not be visible at the leading order iV classical action, but 
only in the order 1 free energy induced from loops. For higher values of h the lowest two 

9 Note that, as the equations of motion only depend on u> 2 , the bands have symmetry around u) = 0. 



14 




FIG. 3: Band structure of the global AdS2 toy model for A = 1 and h = 0.05 (top left panel), 
h = 0.82 (top middle panel), h = 1 (top right panel), h = 5 (bottom left panel) and h = 50 (bottom 
right panel). 

bands undergo an interesting reorganization, changing the topology of the band structure. 
For very large h we once more approach degenerate fc-independent bands (but now with 
u ~ band surviving, as clearly visible in the bottom right panel of Fig. [3]), this time with 
the opposite boundary conditions [in other words, otj has to vanish at each site for h = 
while j5j has to vanish at each site for h = oo, giving t 2 i(cj n ) = rather than f 12 (w n ) = 0]. 10 



V. TOWARDS A HOLOGRAPHIC HUBBARD MODEL 

While the double trace deformations we introduced allow dimer vibrations, or "mesons," 
to propagate through the whole crystal, they do not lead to transport of defect fermions. 
The double trace deformation we have introduced is the product of two U(l)j invariant 
operators, so it preserves all the U(l)j global symmetries. In particular, the difference in 
number of defect fermions at the two ends of each dimer, or "baryon number," is conserved 
and thus the defect fermions cannot move. 

10 Incidentally, the analysis of this section goes through in the same way for fermionic operators Oj. It 
would be interesting to explore implications of this topology change in band structure for such fermionic 
excitations. 
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(a) 

J=l 2 3 




(1,2) (2,3) 

FIG. 4: (a)The high-temperature phase with two branes at each site, with lattice sites equally 
spaced. (b)A possible low-temperature configuration. (c)Another possible low-temperature con- 
figuration on which we will focus. 

The only way to actually introduce a moving charge carrier seems to be to introduce a 
charged field or object in the bulk that can effectively carry baryon number. One potential 
such object would be the FT-boson of a U (2) non-Abelian gauge field living on the stacks of 
branes. So we imagine doing something like: 

1. At each odd J = (2j + l)th site we put two anti-D5-brane, while we put two D5-brane 
at each even J = (2j)th site, and this time we equally space lattice sites (see Fig. 4). 
At each site, we label two species of defect fermions by xl,j and Xr,j- 

2. In order to make charge carriers move around, we consider deforming the Lagrangian 
by a conventional hopping term 

AL hop = t h ^2 x{,jXr,j + Hermitian conjugate. (22) 
j 

Here th can be a complex coupling constant. 

There is a two geometrically distinct configurations with the same energy (see Fig. 4). De- 
formations of the theory as well as l/N corrections can lift this degeneracy. Here we will 



16 



exclusively focus on the state depicted in Fig. 4(c). 

The scalar partner of the W-boson is dual to an operator of the form 



X l ,jXr,j on a given 



site, so turning on a nontrivial source for such a field activates a conventional hopping term in 
the system, rather than the terms quartic in defect fermions we have been implicitly dealing 
with by resorting to double trace deformations. This hopping operator affects the dynamics 
at leading order and is not suppressed like the double trace deformations we considered 
above. Of course, one can always choose th to be a small parameter, and treat the problem 
perturbatively in t h for \t h \ <C 1. 

In the bulk this operator maps to an off-diagonal component of the U (2) j current living on 
the D-branes on the J-th site. Note that in the full D3/D5 system the coupled scalar/vector 
sector on the J-th site is dual to the operators x a jXb,j, the chiral condensate, and x a jloXb,j, 
the dual current. As 70 = i is just a number, these two operators are really just real and 
imaginary parts of the hopping term deformation on each site. The labels a, b run over L 
and R. There are 4 real operators worth of terms we can add to the Lagrangian, and there 
are 4 real operators worth of conserved currents dual to the massless gauge field mode in 
the bulk. As discussed before, in the full D3/D5 system there is a second operator with 
the same global charge assignments and dimension 3. For minimalistic bottom- up models 
without such an operator presumably no bulk scalar field mixing with the vector is required. 

In order to analyze the effect of turning on the off-diagonal components it is convenient 
to treat the problem as a simple U(2) Nhs - gauge theory living on a segment of AdS 2 , where 
iVi. s . denotes the number of lattice sites. Let us parametrize the gauge field on each site as 



A convenient gauge choice is A Z: j = 0. For a minimalistic model one can take the intrinsic 
metric to be a segment of pure AdS 2 



running from the UV boundary at z = to some "hard wall" at z = z*. The hard wall at 
z* simply reflects the fact that we are studying a state of the system in which the branes 
are reconnected, and do not reach all the way to the horizon. The detailed state of the field 
theory is captured by the boundary conditions imposed at the hard wall. 




(23) 




(24) 
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For the purpose of phenomenological model building, the boundary conditions imposed 
in the IR are part of the input. But there is a particular "geometric" set of IR boundary 
conditions that corresponds to the connected bridge configurations we studied before. There, 
the boundary conditions on the diagonal components of the gauge fields follow by continuity 
of the fields and their radial derivatives in the true bridged configuration. If we let a denote a 
single valued radial coordinate along the brane, with z(a) an increasing function on the right 
brane and a decreasing function on the left brane, then the natural boundary conditions are: 

4 L j+i-^5 = °> d a (At + i - = d z A^ +1 + d z A*j = atz = ^. (25) 

Note that these boundary condition gives the desired breaking of the U (2) Nhs - to the U(l) Nhs - 
associated with the U(l) gauge fields living on the bridges: 

= A t = A w+iy ( 26 ) 

Having different boundary conditions on the left and right fields directly breaks each U(2) 
to U{\) x U{1) and then, as in the brane setups studied earlier, the boundary conditions 
ensure that the left U(l) field of the (J + l)th site is identified with the right £7(1) field of 
the Jth site. 

Last but not least, we have to determine the boundary conditions on the off-diagonal 
components of the gauge field, ht,j and h* t 3 . From the bottom-up point of view the best 
way to think about the IR boundary conditions is to introduce an "IR-brane-localized Higgs 
field." As in the cases of interest the boundary conditions always connect U(2) gauge fields 
living on neighboring sites, one can introduce complex, bi-fundamental scalars connecting 
neighboring sites. That is, for every bridge with label (J, J + 1), one adds an IR brane 
localized Lagrangian for a 2 by 2 matrix of scalar fields <f>(j t j+i) 

Tr | 2 = Trldpfaj+i) - iA^j<p {J j +1) + i(j)^j + x)A^j +x \ 2 . (27) 

One particularly interesting form of the vacuum expectation value (vev) of the IR-brane 
localized Higgs field is 





(J.j+i) 



m 



(28) 



Eq.(28) is the unique choice if we want a vev that only gives quadratic terms mixing 



Afj with Af J+1 but no other quadratic terms involving the diagonal components of the 
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gauge fields. This form of the vev is "geometric" in the sense that it can describe branes 
reconnecting as we saw above. This form of a scalar field expectation value gives rise to IR 
brane localized mass terms for the gauge fields of the form 

Lm = £H 2 - Alf + |M 2 + \Kj + i?) ■ (29) 

j 

The boundary conditions on a gauge field with a finite boundary mass matrix (m 2 ) ab (due 
to the boundary Higgs) are in general (32j E3] 

d z A a t = (m 2 ) ab A b t at z = z*. (30) 



We see that the above Higgs vev and the resulting IR Lagrangian (29) in the limit of large 



m give exactly the boundary condition (25) that we know encode the correct geometric 



conditions on the bridges, together with a Dirichlet boundary condition on the off-diagonal 
components of the gauge field. 

In this language it is now straightforward to turn on the actual fermion hopping interac- 



tion of Eq. (22). Asymptotically, we have 

^'j L = af L + (3j> L /z, h tJ = rj + 5j/z. (31) 

Turning on a nontrivial hopping interaction simply tells us that we are studying bulk gauge 
field configurations in which we impose the UV boundary condition that 7 j = th- 
in the full D3/D5 system it is impossible to study this deformation reliably. For one 
thing, turning on the off-diagonal gauge field components requires one to know the full 
non-Abelian DBI action governing the dynamics of these fields. However, this action is not 
known beyond the few lowest dimension terms in powers of F^ u . Even worse, to reliably 
study the bridged configurations we really need the analogue of the DBI action that governs 
a brane/anti-brane system including the tachyon field. This is certainly beyond the scope 
of present-day D-brane technology. 

On the other hand, at the level of bottom-up model building, we have assembled all of 
the ingredients we need to study a holographic realization of a generic lattice model with 
hopping fermions. Thus, one can take one's favorite unsolved lattice model (for example 
the Hubbard model), and parametrize one's ignorance by writing down a corresponding 
higher dimensional brane system with an effective action for both the gauge field and the IR 
Higgs field. This action will have free parameters, which should be matched against known 
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properties of the boundary lattice model. While it is not obvious that such a rewriting 
will be advantageous, it may offer some new approaches to this class of problems, just as 
bottom-up models of hadron physics have done for the study of QCD. 
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In this section we present a simple toy model of probe brane fluctuations in the dimerized 
phase. On the probe brane, let us suppose that there is a scalar field governed by the 
following effective action: 
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Appendix A: Hypergeometric-ology 




(Al) 



where is a caricature "two-AdS2" induced metric given by 



(caricature) ^ ^ 



1 



2 (-dt 2 + dx 2 ) , 



x E [0,a]. 



(A2) 



[cos {% ( 



X 



a 



I)}] 



This describes a global AdS 2 on the bridge. 

Making a coordinate transformation to p = 7r(| — |), we get: 




(A3) 
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By a further transformation 



4>M 



1 + ship 



1 



, «e[0,i], 



v(l — v) 

the equation of motion can be brought into hypergeometric form: 



d 2 ^ 



+ 



dv 2 4v 2 (l-v) 2 



1-4 



au 



Tt 



l-4l™VL+lf-(^)} 



/CLOU 
\ 71 



4>u 



We can now bring this equation into the standard Gauss' hypergeometric form: 



ypcr 



au 

1 

V 



+ r 



''hyper 



au 



7T 



1 /am\ 2 3 

+ 2' C hypcr( c hypcr — 2 J = J ~ ^ 



hyper 

~~2 



c hyper a hyper ^hypcr ^ 
2 



d 2 d 

f(l f J ('-'hyper (^hyper "I - ^hyper "I - l)^) Ohypcr^hypcr 



(A4) 
(A5) 

= 0. 
(A6) 

(A7) 

(A8) 
F u (v). (A9) 



After reorganizing a bit we get 11 



A(A-l) 

<f>(v) 



/ma\ 



y — J [choose the positive root so that (A — 1) > —A] 
1 



l-v 



A-l 
2 



[a(u) 



F 1 1 - A + 



aw 



7T 



. ,1 /aw\ 1 /au;\ 1 

+«-)^^(2 + (t)-2-(t) ; 2 + a -" 



1 



1 - u 



A-l 
2 



1 — Zl\ 

' ^ A + (^),A-(^);i + A;l 



The hypergeometric functions appearing above are connected through frequency- 
dependent matrices: 12 

2A-1 

(3 " ^(l-A+^l-A-^f-A;*;) 
t 12 {u) t 22 {u) \( 2 F, (§ + (f ) , | - (<f ) ; | - A; 1 - „) 

I 1-2A 

t 2 iM M (i) 2 2 F 1 (A+(f),A-(f);I + A;l-t;) 



11 Note that u ~ (p + |) 2 as p — > — | whereas (1 — v) <~ (p — |) 2 as p — > |, giving rise to | in the exponents. 

12 They are also A-dependent, where A is a parameter of the theory. 
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where Uj(u) can be expressed in terms of Gamma functions as 

t u (u) 

*2i(w) 
and t 2 2{oo) 
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